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Abstract 

The hadronic decays B — > J/ipK(K*) are interesting because experimentally they are the only 
color-suppressed modes which have been measured, and theoretically they are calculable by QCD 
factorization even the emitted meson J/ip is heavy. We analyze the decay B — > J/tpK within 
the framework of QCD factorization. We show explicitly the scale and 75-scheme independence of 
decay amplitudes and infrared safety of nonfactorizable corrections at twist-2 order. Leading- twist 
contributions from the light-cone distribution amplitudes (LCDAs) of the mesons are too small 
to accommodate the data; the nonfactorizable corrections to naive factorization are small. We 
study the twist-3 effects due to the kaon and find that the coefficient a,2(J/ipK) is largely enhanced 
by the nonfactorizable spectator interactions arising from the twist-3 kaon LCDA 0^, which are 
formally power-suppressed but chirally, logarithmically and kinematically enhanced. Therefore, 
factorization breaks down at twist-3 order. Higher-twist effects of J ftp are briefly discussed. Our 
result also resolves the long-standing sign ambiguity of 02 (J /ipK), which turns out to be positive 
for its real part. 
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I. INTRODUCTION 



There are several reasons why the decays B — > J/ipK and J/ipK* are of great interest. 
Experimentally, they are the only color suppressed modes in hadronic B decays that have 
been measured so far. These decays receive large nonfactorizable corrections, or equivalently, 
the conventional parameter a 2 (J/ipK) is large, of order 0.20 — 0.30 JT|. In principle, the 
magnitude of analogous a 2 also can be extracted directly from the decays B — > D(*> 7i°(p°) 
and indirectly from the data of B~ — > D^7f(p) and 5° — > D^n(p). However, the former 
color-suppressed decay modes of the neutral B meson are not yet measured. Besides the form 
factors, the extraction of a 2 from B — > D^ir(p) depends on the unknown decay constants 
fo and fr>*. On the contrary, the decay constant fj/^ is well determined and the quality of 
the data for B — > J/ipK^ has been significantly improved over past years. 

From the theoretical point of view, the prominent question is how to calculate the param- 
eter a 2 {J/'4)K). In the literature, this decay mode has been calculated using QCD sum rules 
and the hard scattering approach. Khodjamirian and Riickl have applied light-cone sum 
rules to study the nonfactorizable effects in B — > J/ipK s and concluded that a 2 (J/ipK) is 
negative, whereas Li and Yeh |§ found a positive a 2 (J/ipK) based on the perturbative QCD 
hard scattering approach. Therefore, there is a sign ambiguity for a 2 ( J/ipK). Although it 
has been argued that a negative a 2 (J ' ji\)K) is very unlikely for several reasons |TJ], it is impor- 
tant to have an independent theory calculation to clarify the sign issue. The QCD-improved 
factorization approach advocated recently by Beneke, Buchalla, Neubert and Sachrajda j| 
is suitable for this purpose. In this approach, nonfactorizable effects in B — > M\M 2 with 
recoiled M\ and emitted light meson M 2 are calculable since only hard interactions between 
the (BMi) system and M 2 survive in the heavy quark limit. 

The aforementioned QCD factorization method is no longer applicable if the emitted 
meson is heavy. For example, since the D° meson produced in B° — > tt°D° decay is not a 
compact object with small transverse extension, it will interact with the (Btt) system in the 
presence of soft interactions. In principle, the parameter a 2 (irD) cannot be calculated using 
QCD factorization, though it has been roughly estimated in M by treating the charmed 
meson as a light meson, which is certainly a dubious approximation. Fortunately, the QCD 
factorization approach can be applied to B — > J/tpK decay since the transverse size of J/ip 
becomes small in the heavy quark limit. However, a recent study by Chay and Kim || 
indicated that while the factorization method is applicable to B — > J/ipK, the leading-twist 
contributions are too small to explain the data. 

In general, power corrections to QCD factorization are suppressed by a factor of Aqcd/^16- 
However, as shown in || and Q, there exist power corrections which are chirally enhanced. 
Moreover, there are also some power-suppressed terms at twist-3 order which involve the 
logarithmically infrared divergence, indicating the dominance of soft gluon exchange. That 
is, there are twist-3 power corrections which are chirally and logarithmically enhanced. In 
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the present paper, we find that a twist-3 light-cone distribution amplitude of the kaon can 
lead to a large spectator interaction, which allows to alleviate the discrepancy between theory 
and experiment for B — > J/ipK. 

The purpose of this work on B — > J/ipK within the framework of QCD factorization is 
twofold. First of all, we perform similar leading- twist analysis as in || . Second, we investigate 
higher- twist effects on a 2 (J/if)K) and show that factorization breaks down. The rest of this 
paper is organized as follows. In Sec. II, we review the QCD factorization approach and 
introduce twist-2 light-cone distribution amplitudes of mesons. Then we proceed to compute 
vertex and spectator corrections to B — > J/ipK. In Sec. Ill we discuss the corrections at 
twist-3 order and give numerical results and discussions in Sec. IV. The conclusion is given 
in Sec. V. 



II. B -> J/ipK AT LEADING-TWIST ORDER 

A. QCD factorization 

Consider the hadronic decay B — ► M\M 2 with Mi being recoiled and M 2 , which is a 
light meson or a quarkonium, being emitted, it has been shown that the transition matrix 
element of an operator O in the effective weak Hamiltonian valid up to corrections of order 
^QCD/ m b is schematically given by || 



{M 1 M 2 \O l \B) = (M 1 \j 1 \B){M 2 \j 2 \0) 



l + 5> n a? + 0( 



2 Ff M \m\) jf 1 duT^u^MM 



^QCD- 



+ / rf^M^T // (^,M,u)0 i? (^)0 Ml (f)0 M2 (M), (2.1) 







where F BMl is a B — Mi transition form factor, <pM is the light-cone distribution amplitude, 
and T , T n are perturbatively calculable hard scattering kernels. In the naive factorization 
approach, T 1 is independent of u as it is nothing but the meson decay constant. However, 
large momentum transfer to M 2 conveyed by hard gluon exchange implies a nontrivial convo- 
lution with the distribution amplitude 4>m 2 - The second hard scattering function T n , which 
describes hard spectator interactions, survives in the heavy quark limit when both M\ and 
M 2 are light or when M\ is light and M 2 is a quarkonium ||]. 

The factorization formula (|2.1| ) implies that naive factorization is recovered in the 
rrib — > 00 limit and in the absence of QCD corrections. The radiative corrections to naive 
factorization are calculable since soft gluon interactions between the {BM\) system and the 
M 2 meson are power-suppressed in the heavy quark limit. Consequently, the nonfactorizable 
contributions to naive factorization are actually amenable in the infinite quark mass limit. 
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The QCD factorization formula is not applicable to the decay, for example, B° — > tt D 
where the emitted meson D° is heavy so that it is neither small ( with size of order 1/Aqcd) 
nor fast and cannot be decoupled from the {Bit) system. This is also ascribed to the fact 
the soft interaction between (Bit) and the c quark of the D° meson is not compensated by 
that between (Bit) and the light spectator quark of the charmed meson. 

In Sec. II. C we show that the QCD factorization formula is still applicable to B — > J/if)K 
decay though the emitted J / if) is heavy. The point is that the transverse size of J / if) becomes 
small [of order l/(m c a s )} in the heavy quark limit. Technically, infrared divergences arising 
from the soft interactions between the c quark of J /if) and {BK) system and between the c 
quark and {BK) compensate. Therefore, the nonfactorizable contributions to B — > J/if)K 
is infrared safe. 

In principle, power corrections are of order Aqcd/^6 and hence they can be neglected 
in the heavy quark limit. Nevertheless, there are some power corrections which can be 
enormously enhanced and hence cannot be neglected. First of all, the contributions to, for 
example B — > Kit, from the (5* — P)(S + P) penguin operators are enhanced by the factor 



2fi x _ 2m 2 K _-Mqq) in A 



12 ^£ (2.2) 
m b (m s + m u )m b m b ffc m b 

which is proportional to the quark condensate. Second, it is shown in || that the hard 
spectator interaction in B — > Kn to twist-3 order has the form 



, 4tt 2 f K f B r 1 dp B r 1 d£ K r 1 dr] r 2 // x 

flI = ir c WW^ B L {p) l 1^ {0 l j r {v) + ^ {v) 



(2.3) 



Since the twist-3 distribution amplitude <pp(fj) ~ 1, it does not vanish at the endpoints. 
Consequently, the logarithmic divergence of the fj integral implies that /// is dominated by 
soft gluon exchange between the spectator quark and quarks that form the emitted pion, 
indicating that factorization breaks down at twist-3 order. Hence, this power correction 
is chirally and logarithmically enhanced. We shall see in Sec. Ill that the hard spectator 
interaction in B — > J/if)K receives the same chirally enhanced infrared logarithms from 
twist-3 kaon distribution amplitudes. 

B. Light-cone distribution amplitudes of mesons 

Consider the matrix element of nonlocal operators sandwiched between the vacuum and 
the vector meson J/if>: 

(J/iPK(x)cl(0)\0) = |l{(J/^|c(x)c(0)|0) + 75 (W|c(x)75c(0)|0) +7^(W|c(a;)7 M c(0)|0) 



4iV c 

-7^75 W|c(x) 7/ ,75C(0)|0> + \(T^{J/il>\c(x)(T^c{0)\0)} (2.4) 
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where a, b are color indices, a, (3 are indices for Dirac matrices. The leading-twist light-cone 
distribution amplitudes (LCDAs) of J/ip are given by M 

C*W . X fl 111 

(W(P,A)|c(x) 7 ,c(0)|0) = /^m^-^—P^jf e^"^(0, 

(J/^(P, A)|c(x)^c(0)|0) = -ifJ/^P, -e^P^J^e^^), (2-5) 

where e* is the polarization vector of J/ip, £ is the light-cone momentum fraction of the c 
quark in J/i/j, fj/^p and fj,^ are vector and tensor decay constants, respectively, but the 
latter is scale dependent. The normalization conditions of the twist-2 LCDAs are 

/ 1 ^0^(O= / X de^(0 = l- ( 2 -6) 
Jo 11 jo 

Likewise, at the twist-2 accuracy, the only kaon LCDA which contributes to the matrix 
element (the relation o^ v <g> a^ v = a^ u -y 5 ® cr^vlb being applied) 

Sjab 

(K-\s a a (0)4(x)\0) = — {<JT-|5(0)«(x)|0> + l5 (K~\s(0) l5 u(x)\0) +y*(Jir|a(0)7 M «(x)|0> 

-7^ 75 (^|s(0) 7m 75w(x)|0) + ^^(K-I^O)^^^)^)}^, (2.7) 

is given by |§ 

(A-(P)|S(0) 7M 7fe«(a;)|0> = -i/^ f dfje^^if)), (2.8) 

Jo 

where 77 is the momentum fraction of the antiquark u in K~ and Sl(j) K {f])dfi = 1. 
As for the P meson LCDA, we will follow [|J to choose 

(0|g Q (x)6^(0)|75(p))U +=x±=0 = -!^[(^ + mB ) 75 ]^^^pe-^-[0f(p) (p)] 7Q , (2.9) 

based on the observation that the B meson is described by two scalar wave functions at the 
leading order in l/m&, where p is the momentum fraction carried by the spectator quark of 
the B meson. In Eq. ( |2.9|) , n_ = (1,0, 0,-1) and the normalization conditions are 

f 1 dptf(p) = 1, f dptfifi) = 0. (2.10) 

Jo Jo 

The leading-twist LCDAs of J/ip can be expanded as M, 

0^(0 = sea - (1 + \ 4 im - 1) 2 - 1]) , 

= 6£(1 - (l + \ 4 [5(2£ - l) 2 - 1]) , (2.11) 
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where the parameters a| and are defined by the matrix element of a twist-2 conformal 
operator with conformal spin 3 [TAP] while twist-2 DA (p K can be expanded in terms of 
Gegenbauer polynomials CfJ 2 |§: 



in^ 2 ) = 677(1 -fj) ( 1 + y: aU^clL 2 ^ - 1 



(2.12) 



n=l 



with the values of the Gegenbauer moments a% being available from ||. For the B meson, 
we use [p] 



4>f(p) = N B p 2 (l-p) 2 exp 



1 / pm B 

2 V oj b 



(2.13) 



with uj b = 0.25 GeV and N B being a normalization constant. This B meson wave function 
corresponds to X B = 303 MeV defined by Jq dp 4> B (p) / p = m B / X B . This can be understood 
since the B meson wave function is peaked at small p: It is of order m^/AgcD at p ~ 
Aqcd/^b- Hence, the integral over 4> B (p)/p produces an m B /Aq CD term. 



In ensuing calculations we find that the contributions from the J/ip LCDA 



,Jli> 



arc 



proportional to the factor (/jy^^ c )/(/j/»/) m J/»/>) relative to that from 0j[ • Contracting 
(J/'0|c(O)a' Ati/ c(O)|O) in Eq. (2.4) with P u and applying the equation of motion and the 
matrix element 



we are led to 



Therefore, 



(J/V|c(0)7 M c(0)|0) = fj/^mj/^e*, 



(2.14) 



(2.15) 



where in the last step we have applied the on-shell relation £ Fjm 
(not anticharmed) quark in the J/ip- 



(2.16) 

rrii for the charmed 



Since the parameters a| and are unknown, we will employ the asymptotic DAs </>j[ (£) 



1 

Sec. IV. For general discussions in Sec 

III, we still use Eq. (PU) for the LCDAs of J/ip. 



±^{C) = 6£(1 — £) in numerical calculations in Sec. IV. For general discussions in Sees. II and 
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C. B J/^lf in QCD factorization 



The effective Hamiltonian relevant for B — > J/ipK has the form 



7i 



cff 



where 
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cMOM + c 2 (/i)o 2 (/i) - v^; £ Ci(n)Oi(n) + h.c, (2.17) 



i=3 



Ci = (cfo) v _ A (sc) v _ A , 

03(5) = (^l^E^V-il 



(4 



(v+a)j 



7(9) 



( s6 )v-AE e 9'(?V) 



04(6) = (Sabp) v _ A J2(<l , p<l'a)v-A(v+A)i 
q' 

3 



o, 



8(10) 



(2.18) 



v+a(v-a)? 



q> 



q' 



with 3 -0 6 being the QCD penguin operators, 7 -O w the electroweak penguin operators, 
and {q x q 2 ) = qa^(l - 7s)<?2- 

Under naive factorization, the decay amplitude of B — > J/ipK reads 



Gi 



A(B -> J/Vlf) = -^V c6 V;;(a 2 + a 3 + a 5 + a 7 + ag^mj^Ff K (m^)(2e* ■ p B ), (2.19) 

where a 2 j = c 2 j + (l/AQc^-i and a 2 i-i = c 2 j_i + (1/-/V" c )c 2 j in naive factorization and the 
approximation V t bV t * s ~ —V c bV* s has been made. The form factor is defined by 



rrip — m P , 



mp — mp 



Q 



-q,F Q (q 2 ). (2.20) 



There are two serious problems with the naive factorization approximation. First, the Wilson 
coefficients and hence a, are renormalization scale and 75-scheme dependent, whereas 
the decay constants and form factors are not. Hence, the amplitude ( |2.19 ) is not physical. 
Second, nonfactorizable effects, which play an essential role in color-suppressed modes, are 
not taken into account. 

The aforementioned two difficulties for naive factorization are resolved in the QCD factor- 
ization approach in which the inclusion of vertex corrections and hard spectator interactions 
(see Fig. 1) yields^ 



ci a s C F 
a-2 = c 2 + — + — — ci 

N n 47T N r . 



18 
14 



19 In ^ 4- f 4- f 4 I !*0<™mL n 
iZ[n r JI + Jll + t-iRK / 3 \9l 



^Using the constant matrices tndr and fnv given in Eqs. (2.18) and (2.19) in |l0|, it is straight- 
forward to obtain the constant terms in Eq. ( |2.21 ) in NDR and HV 75-schemes. 
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a, CV 



a3 = C3+ iV^ + 4^iV: C4 



12 In 



a 5 



c 6 



C5 + A^-4^i^ C6 



JL 

m b 



fi + fn + 



F 



BK i 



rn 



Jli» 



J/i> 



9i 



o 



121n-^ + / J + / /I + 
nib r l 



a 7 = c 7 + 



C8 



a; ., Cf 



JV r 4tt JV, 



■c fl 



, Cio , a s C F 



i2i n A + /j + /w + 

- 121n-^- + // + /// + 
m 6 



J/i" 



Ff K (m 



9i 



9i 



(2.21) 



tiBK, 



111 



2 



F/^(m 2 



where the upper entry of the matrix is evaluated in the naive dimension regularization 
(NDR) scheme and the lower entry in the 't Hooft-Veltman (HV) renormalization scheme, 
Cp = (N^ — 1)/(2N C ), and N c is the number of colors. The hard scattering functions // 
and gi arise from the vertex corrections, Figs. l(a)-l(d), while /// from the hard spectator 
interactions Figs. 1(e)- 1(f). Formally, the coefficients are scale and 75-scheme independent 
and we will come to this point again in Sec. IV. 

The results for the hard scattering functions // and gi are 



Si 




(2.22) 



and 



9i 



4£(2£-l) 



+ 1 



1 



z)(l 
1 







ln£ + 



[1 



z(l 



ln(l - z) 



1 - zO 2 [1 
2(l + z-2z£) 



z£ In z£ 



VK 



Z)(l 



[l-z(l -0] s 



(2.23) 



where z 



m 



J/i> 



/m 2 B and the assumption J// ^(£) = <fif/^(£) = <j>±^(0 has been made. In 



deriving Eqs. ( gjp and ( ggg) we have applied (|T|) and the relation J/ ^(O = J/V> (O, 
where £ = 1 — £ is the light-cone momentum fraction of the c quark in J/i/j. 

Several remarks are in order, (i) We have proved explicitly the cancellation of infrared 
divergences so that the resultant amplitude is infrared finite when summing all the diagrams 
in Fig. 1, a key element for the applicability of QCD factorization. Moreover, we also 
show a cancellation of the infrared double poles, i.e. l/e 2 R , as in pT| . As stressed and 
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c c c c 




to (/) 



FIG. 1. Vertex and spectator corrections to B — > J/ipK. 

elucidated in H, in the limits 2: ^ and z / 0, the amplitudes are infrared finite in 
both cases, but for different reasons, (ii) The strong phases in // and gi arise from the 
diagrams Figs. 1(c)- 1(d) where there are hard gluon exchanges between the outgoing K 
and J/ip. (iii) Only the form factor Ff K contributes to the decay amplitude under naive 
factorization, while to the order a s both form factors Fq and Ff K contribute, (iv) Our 
results for // and gi agree with that in || by noting that the form factors Fq± are related 
via Fq K (m 2 ,^) / Ff K (m 2 j^) = {m 2 B — m 2 ^) / m 2 B p. The only difference is that we treat the 
ratio (fJ/^Tn c ) / (fj/^mj/^) as 2£ 2 [see Eq. (|2.16|) 1, while it is considered to be a constant in 
||. (v) It is easily seen that in the zero charmed quark mass limit, 

fj + gj = f* W+it) ^i=^ln£-3i7r) , (2.24) 

in agreement with for B — > tttt, as it should be. 

As for the hard scattering function /// originating from spectator diagrams, we write 

fn = fh + f!i + ---, (2-25) 

where the superscript denotes the twist dimension of LCDA. To the leading-twist order, we 
obtain 
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47T 2 /y/ 

iV c F 1 BK (m 2 J/ ^ml Jo 

x ^i+(p-2e+^+4^ (2>26) 

p(p-v + vz)[(p ~0(p-v) + (vp pO z \ 

This can be further simplified by noting that p ~ C(AqcD/ m f>) —> in the — > oo limit. 
Hence, 



/? , = £ fefe , j_ fdptM. f M tarn f in m, (2.27) 

A/ c Ff A {mj^)m B 1 — zJo p Jo £ Jo r\ 

where the z terms in the numerator cancel after the integration over 

Before proceeding we make two remarks: (i) At a first glance, it appears that the spectator 
contribution is power suppressed by order (Aq CD /m b ) 2 . This is not the case. First of all, the 
B wave function <pf(p) is of order m&/AQCD at p ~ Aqcd/^& [see also the discussion after 
Eq. (p.9|)1- Second, the form factor Ff K scales as ~ (AQCD/^b) 3//2 , while decay constants as 
Jb ~ Aqcd / m t an< ^ /?r ~ Aqcd 0. It follows that ffj is of the same order as // and gi and 
is not power suppressed, contrary to the claim in || . Furthermore, the former is numerically 
more important than the latter, as we shall see in Sec. IV. (ii) By inspection of fl2.26|) , there 
will be some strong phases coming from hard spectator interactions beyond the heavy quark 
limit, which are difficult to evaluate numerically, however. 



III. HIGHER-TWIST EFFECTS 



In the last section we have computed leading-twist vertex and hard spectator corrections 
to naive factorization. However, we shall see in Sec. IV that numerically the twist-2 nonfac- 
torizable effects are small; the predicted decay rate of B — > J/tpK is too small by a factor 
of 7 ~ 10. Therefore, it is inevitable that higher-twist effects which are seemingly power 
suppressed should play an essential role. 

As shown in Sec. II. A, twist-3 power corrections to the spectator interactions could be 
enormously enhanced by the chirally enhanced infrared logarithms. One reason that the 
spectator diagrams could receive large power corrections is ascribed to the fact that the hard 
gluon exchange in the spectator diagram is not hard enough. The virtual gluon's momentum 
squared is 

k 2 = (-pp B + Vp K ) 2 ~ ~PV m B ~ -A QC D^b, (3.1) 



^"Although the same expression ( 2.27 ) is also given in ||, it is not clear to us how to achieve 
the cancellation of the z terms appearing in the numerator of Eq. ( 2.26Q when 4£ 2 is replaced by 
2(fJ/- t p m c)/(fj/ip'mj/ii,) and treated as a constant according to ||. 
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which is of order 1 GeV due to the smallness of the momentum fraction p ~ AqcdA 7 ^ 
carried by the spectator quark in the B meson. Therefore, we shall study higher-twist power 
corrections to hard spectator interactions. 

A. Twist-3 LCDAs of the kaon 

We first consider the twist-3 DAs <ft K and (ft K of the kaon defined in the pseudoscalar and 
tensor matrix elements ||: 

(lT(P)|s(0)<7bt»(a:)|0> = fK f K /* dfje^fifj), 

m s + m u Jo 



6 m s + m u 



m s + m u x - 
m K 



X(P^ U - P u x,) f dfje if i p - x ^(fj). (3.2) 

J 



They can be expanded in terms of Gegenbauer polynominals 

^(r 1 ) = l + aC 1 2 /2 (r ] ) + bCl /2 (r ] ) + ... 



v 

tfirj) = 6n(i -m + dcTin) + ■■■), (3-3) 



where the coefficients a, b, d can be found in j8|. From Eq. ( |3.2|) it is clear that twist-3 
DAs of pseudoscalar mesons are associated with a chiral enhancement factor p x . As before, 
fj is the light-cone momentum fraction of the u quark in K~. Just as the hard spectator 
interactions in B — > Kir decay which receive large power-suppressed corrections from twist-3 
LCDAs of the pion (see Sec. II. A), we find that the twist-3 kaon LCDA <ft K contributes to 
spectator diagrams in B — > J/iftK decay and yields 

l\ c t 1 {rrij^)m B [m s + m u )m B Jo 

x 1 [ 1 AZ + 8 ^ \ (3 4) 

6 \£pfj 2 (l — z) 3 pfj 3 (l — z) 3 pfj 3 (l — z) 3 J 

Applying the twist-3 DA <ft K {ff) — 677(1 — fj), we see that the linear divergent terms propor- 
tional to Jo dfj/fj 2 cancel after the integration over £; this cancellation happens because the 
J/ijj distribution amplitude is even in (2£ — 1) [see Eq. ( |2.11| )]. Consequently, 



_ /2^\ 47^ f K f B ,1 dp B /•! di j/ip /-i dq (f> K (fj) f3 5) 

}lI -\m B ) N c F?X(m* JM )m* B h p^ {P) Jo £ { ^ } Jo ^6(l-z)»' 

with p x being defined in Eq. (|2.2|) . Therefore, to the twist-3 order of kaon DA, 
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[mj^jmg 1 - z Jo p io ( 




!k!b 1 f 1 dp B f 1 d£ 



h fj \ T v " rriB (1 — z) 2 6fj 



x ^K^ + I^TT^^ > (3-6) 



which, apart from the constant term, agrees with Eq. (|2.3| ) in the z — > limits The 
logarithmic divergence of the fj integral implies that the spectator interaction is dominated 
by soft gluon exchanges between the spectator quark and the charmed or anti-charmed 
quark of J/if>. Hence, factorization breaks down to twist-3 order. Following |J, we treat the 
divergent integral as an unknown parameter and write 

i,r^h(^ +r , (3.7) 



V V^QCD, 

with r being a complex random number. Therefore, although ffj is formally power sup- 
pressed in the heavy quark limit, it is chirally enhanced by a factor of (2/x x /Aqcd) ~ 12, 
logarithmically enhanced by the infrared logarithms and kinematically enhanced by a factor 
of 1/(1 — z) 2 . Numerically, nonfactorizable corrections to naive factorization are dominated 
by ffj (see Sec. IV). 

B. Higher-twist LCDAs of J/ip 

Intuitively it is expected that light mesons produced in hadronic B decays are appro- 
priately described by LCDAs. Since J /if) is not light enough, it is nature to conjecture the 
importance of higher-twist effects of J /if). The twist-3 LCDAs of J /if) are given by [|7| 

(,.*(A) ™ \ „i 
^ (A) "~P^T^) Jo d t^ P "!>* ®> 

1 / 9mA 

(J/if)(P, A)|c(x)7 M 7 5 c(0)|0) = --m J/% b fj/^ - f JN> 



§Just as B — > J/ipK decay, we find that the twist-3 pion distribution amplitude that contributes 
directly to the spectator interactions of B — > Kir is </>J rather than $L. More specifically, we have 
a contribution like X = ^{dfj/fj 2 )^/^. At a first glance, it seems to be different than the last 
term appearing in Eq. (^). However, since the LCDAs 4> p and <f> a are related by equations of 
motion @, <f> ff (fj)/6 = r?(0 p (r?)/2 + ^(r/)/12), it follows that X = J \dfj/fj)(Mn)/^ + VM/ 1 *)- 
This means that X can be expressed in terms of (f> p and cp' a which do not vanish at endpoints. 
Since (f>' a {ff) = 6(1 — 2fj), X is reduced to f Q {dfj / fj)<f)p(fj) after the non-logarithmic term is dropped. 
Nevertheless, we wish to emphasize that the logarithmic divergence arises originally from 4> a instead 

of (j) p . 
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xe^e^PV f dte* p -*gi\£), (3.8) 

J 

for chiral-even DAs g\_ ] and gf' 



(J/V(P, A)|c(a;)^ C (0)|0) = ■ • • - ^(P^ - P^-^—^m 2 ^ j Q d^^h^O, 

(J/V(P, A)|cOr)c(0)|0> = i (fa - f m (e* ■ x)m^ f d£ e« p -fc« (£), (3.9) 

for chiral-odd DAs n» and /i| , where ellipses denote twist- 2 LCDAs given by fl2.5| ). 

The calculation of twist-3 effects for J/ ip is more complicated since it involves four more 

LCDAs and some other technical problems. Nevertheless, we can argue that the DAs g± 

» 



and h<< are less important. From Eq. ( |2.16| ) it follows that the term 



fj/ip ~ ft/* = fj/ ^ ( 1 ~ 5^ i ( 3 ■ 1 ° > 

vanishes in the heavy quark limit. Therefore, it is of order Aqcd/^6- A full study of the 
twist-3 effects of J/ip is in progress. However, just as the B — > Kti decay discussed in Sec. 
II. A, it is expected that twist-3 corrections which manifest predominately in hard spectator 
interactions are dominated by the kaon. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

Before proceeding to the numerical results, we first demonstrate the scale and scheme 



independence of the coefficients a, given in Eq. ( 2.21 ). To show the scale independence of 
03, for example, we note 



din /i ~ 4rr lij 



T^QGu), (4.1) 



where the anomalous dimension matrix 7 can be found in the literature (see e.g. |13|). We 
find 

"( c 3 + t) = 12 — — c 4 (4.2) 



dln/i 3 4vr N c 

and hence das/din^ = to 0(a 2 ), as it should be.Q As for the 75-scheme independence 
of a,i, we take the NLO Wilson coefficients from [T3[] (see Table 1). It is ready to check the 
scheme dependence of a; from Eq. ( |2.21| ). 



**Empirically one will find that there is a slight scale dependence of a, owing to the fact that not 
all leading logarithmic corrections to dj to all orders in a s are included in Eq. (^2l|) . 
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TABLE I. Lowest-order (LO) and next-to-leading-order Wilson coefficients in NDR and HV 
75-schemes at n = mb(mb) = 4.40 GeV for = 225 MeV taken from Table XXII of [13], where 
a is the fine-structure constant. 





Cl 


C2 


C3 


c 4 


C5 


C6 


c-j/a 


c 8 /a 


c 9 /a 


cio/a 


LO 


1.144 


-0.308 


0.014 


-0.030 


0.009 


-0.038 


0.045 


0.048 


-1.280 


0.328 


NDR 


1.082 


-0.185 


0.014 


-0.035 


0.009 


-0.041 


-0.002 


0.054 


-1.292 


0.263 


HV 


1.105 


-0.228 


0.013 


-0.029 


0.009 


-0.033 


0.005 


0.060 


-1.283 


0.266 



From Eq. ( |2.19| ) we see that the B — > J/ipK amplitude is governed by the parameter 

a 2 = a 2 + a 3 + a 5 + a 7 + a 9 . (4.3) 

In naive factorization, it is predicted to be aJj (J/ipK) = 0.074 using the leading-order 
Wilson coefficients at jj, — rrib given in Table I. Evidently it is too small compared to the 
experimental value 

\a 2 (J/^K)\ cxpt = 0.26 ±0.02, (4.4) 



extracted from the data [14] 



B(B~ -> J/ipK~) = (10.0 ± 1.0) x 10" 4 , B(B° -> J/ipK°) = (8.9 ± 1.2) x 10~ 4 . (4.5) 
In ensuing calculations we will use the form factors 

= 0-70, F BK (m^) = 0.50, (4.6) 

from 0, the decay constants f K = 0.16 GeV, f B = 0.19 GeV, A QCD = 300 MeV, and the 



running quark masses [13[]: m;, = 4.40 GeV and m c = 1.30 GeV. For the parameter r in 
( p.7|) , in principle it may be complex due to soft rescattering. In 0, r is chosen randomly 
inside a circle in the complex plane of radius 3 ("realistic") or 6 ("conservative"). 

In deriving the hard scattering functions fj, gi, ffj, ff r , we have neglected the difference 
between m& and in the heavy quark limit. However, since ffj is sensitive to the value 
of z, we will vary z from one extreme z = m^/mg to the other extreme z = m 2 j,,/m 



as 



shown in Table II, where we have used r = 4.5 for illustration. The sensitivity of ffj on z 
is not unexpected because the twist-3 effects on spectator interactions are governed by soft 
gluon exchange. In Table II we show the numerical values of the parameters a2{J/ipK) and 
a,2(J/ipK), and, for comparison, the twist-2 results of a 2 (J/ipK) denoted by a$(J/i/jK), i.e. 
the predictions without the twist-3 effect ffj. 

Since a priori the shape of the J /if} LCDA is unknown, it is worth considering other 
possibilities besides the asymptotic form which we have used thus far. The delta-function 
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form J '^(£) = 8(£ — \)i as suggested in ||, appeals to the naive expectation of the wave 
function. The results are shown in the last row in Table II. We see that while // and gi are 
not sensitive to the change of the J/ip wave function, the values of fh and fh are reduced 
by a factor of 2/3 owing to the fact that the integral fo(d£/£)(p J ^(£) [see Eq. ( |3.6|) 1 is 
equal to 3 for the asymptotic form and 2 for the delta form of the J / ip wave function. As a 
consequence, the resultant \a 2 (J /ipK)\ becomes smaller. 

TABLE II. Numerical results of hard scattering functions and the parameters a^(J/ipK), 
^(J/ipK), a 2 (J/ipK) for two choices of z and for r = 4.5. For comparison, the predictions using 
the delta-function form for the J/if) distribution amplitude are shown in the last row. 





fl 


9i 


fh 


fh 


af(J/^K) 


a 2 (J/^K) 


a 2 (J/i>K) 


Z = m2 J/^/ m B 


-1.17-i6.14 


0.54-i0.74 


4.89 


12.71 


0.051-i0.057 


0.166-i0.056 


0.158-i0.057 


Z = m2 J/^/ m b 


-1.08-i4.60 


0.73-U.27 


6.36 


21.47 


0.065-i0.047 


0.254-i0.046 


0.247-i0.047 


Z = m J/if,/ m B 


-1.03-i6.43 


0.77-i0.79 


3.26 


8.47 


0.040-i0.059 


0.119-i0.059 


0.111-i0.059 



From Table II we see that at twist-2 order, \a^(J/ipK)\ is of order 0.07 ~ 0.08 which 
is very close to the naive prediction a\°(J/ipK) = 0.074. This means that nonfactorizable 
corrections in the heavy quark limit to naive factorization is small. Therefore, the predicted 
branching ratio of B — > J/ipK to the leading-twist order is too small by a factor of 7 ~ 10. 
When the unknown parameter r = \r \ exp(i8) is varied from 3 to 6 in the complex plane and 
z varies from rttji^jrv? B to raj^/mf, we find that \a2(J/4>K)\ and \a2(J/4>K)\ fall into the 
range 0.12 ~ 0.29 for 5 < 45° (see Fig. 2), which is a reasonable choice since the phase of 
the spectator diagrams is expected to be small. 

Needless to say, the unknown parameter X, the dependence of fh on z and the unknown 
J/ip LCDAs are the main sources of theoretical uncertainties. Hence, our result for a 2 {J/ ipK) 
should be regarded as an order of magnitude estimate. Nevertheless, it is clear that the 
discrepancy between theory and experiment is greatly improved by the inclusion of kaon 
twist-3 effects and the sign of Rea2 {J/fp K) is positive. It remains to study the higher-twist 
effect of J ftp. 

V. CONCLUSIONS 

We have studied the decay B —>■ J/ ipK within the framework of the QCD factorization 
approach. Physically, the factorization method is still applicable because the transverse size 
of J/ip becomes small in the heavy quark limit so that its overlap with the {BK) system 
is small. Technically, the infrared divergence due to soft gluon exchange between the c 
quark and the (BK) system is canceled by that between the c quark and (BK) so that 
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0.4 




FIG. 2. The coefficient \a,2(J/i/jK)\ vs. the phase of the parameter r. Solid and dashed curves 
are for \r\ = 6 and 3, respectively. The upper and lower solid curves are for z = fnj^/fnl an d 
m J/i/>/ m B' respectively, and likewise for the dashed curves. 

nonfactorizable corrections are infrared safe. We have shown explicitly that this is indeed 
the case to the twist-2 order. 

The scale and 75-scheme problems with naive factorization are resolved when radiative 
vertex corrections to the hadronic matrix elements are included, as we have proved explicitly. 
However, nonfactorizable corrections to naive factorization due to vertex diagrams and hard 
spectator interactions are small. Therefore, the predicted branching ratio of B — > J/ipK to 
the leading-twist order is too small by a factor of 7 ~ 10 compared to experiment. 

Since the virtuality of the hard exchanged gluon in the spectator diagrams is only of order 
?™&Aqcd, it is expected that power- suppressed terms arising from higher-twist wave functions 
become important. We found that the contribution from the twist-3 kaon light-cone distribu- 
tion amplitude ^ to spectator interactions is not only chirally but also logarithmically and 
kinematically enhanced, indicating that soft gluon exchange between the spectator quark of 
the B meson and the charmed or anticharmed quark of the charmonium dominates. Conse- 
quently, factorization breaks down at twist-3 level. In spite of many theoretical uncertainties, 
the predicted 1 02 {J/ip K)\ can range from 0.12 to 0.29. This also solves the long-standing 
sign ambiguity for Kea2(J/ipK). Though it remains to check the higher- twist effects of the 
charmonium, it is expected that twist-3 corrections which manifest predominately in hard 



16 



spectator interactions are dominated by the kaon. 
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